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Abstract. We consider a system of N bosons interacting through a two-body potential with, pos- 
sibly, Coulomb- type singularities. We show that the difference between the many-body Schrodinger 
evolution in the mean-field regime and the effective nonlinear Hartree dynamics is at most of the 
order 1/A'^, for any fixed time. The Af-dependence of the bound is optimal. 

1. Introduction 

A system of N bosons in three dimensions is described by ipN G -^s(^^^! 'i^i ■ ■ ■ ^xn), the subspace 
of L'^{M.^^ , dxi . . .dx^) consisting of functions which are invariant with respect to permutations of 
the N particles (the invariance w.r.t. permutation expresses the bosonic symmetry; fermionic systems 
are described by antisymmetric wave functions). We always assume ip^ to be normalized so that 
llV'AfIb — 1 {\'4'n{xi, . . . ,xn)\'^ is interpreted as the probability density for finding particles close to 
(a;i, . . . ,xjv))- We consider Hamilton operators with two-body interactions, having the form 

N N 

hn,x^J2-^-^ + ^J2^^^^-'^o) (1-1) 

and acting as self-adjoint operators on the Hilbert space L'^{M.^^ , dxi . . .dxN)- In the sum of 

the Laplacians is the kinetic energy of the N particles, A G M is a coupling constant, and the sum of 
V{xi — Xj) over all pairs of particles describes the potential energy of the system {V{xi — Xj) acts as 
a multiplication operator); appropriate conditions on V will be specified below. 
The evolution of the system is governed by the N particle Schrodinger equation 

idtil^N.t = HN,\il^Nd ■ (1-2) 

The solution of the Schrodinger equation can be obtained by applying the unitary group generated by 
Hn,\ to the initial wave function ipN,t=o', in other words, (jl.2p is always solved by iJjN,t = e~'^™'^*'0Ar,o- 
In this sense, establishing existence and uniqueness of solutions of (jl.2[) is not an issue. What makes 
the study of (|1.2p challenging is the fact that, in systems of interest in physics, the number of particles 
N involved in the evolution is typically huge {N ranges from values of the order 10'^ in extremely 
dilute samples of Bose- Einstein condensates, up to values of the order 10^'^ in chemical samples). For 
such values of N, the expression ipN,t = e~'^"'^*V'Af,o is not useful if one is interested in establishing 
quantitative or even qualitative properties of the dynamics. For this reason, one of the main goals of 
quantum statistical mechanics is the derivation of effective evolution equations which, on the one hand, 
can be approached by numerical methods (in contrast with (jl.2l) ). and, on the other hand, approximate 
the solution of (|1.2p in the interesting regimes. 

One of the simplest regime where effective evolution equations can be used to approximate the full 
many-body evolution is the so-called mean field limit, which is characterized by a large number of 
very weak collisions among the particles. To realize the mean field limit, we consider large values of 
N (many collisions) and small values of the coupling constant A (weak interactions). A non-trivial 



Li Chen is partially supported by National Natural Science Foundation of China (NSFC), grant number 10871112. 
Benjamin Schlein is partially supported by an ERC Starting Grant. 

1 



2 



LI CHEN, JI OON LEE, AND BENJAMIN SCHLEIN 



effective dynamics can only emerge when the many collisions produce a total force of order one on each 
particle; in other words, when NX is of order one. To study the mean-field regime, we set therefore 
A = 1 /N and we consider the evolution generated by the Hamiltonian 

N ^ N 

i/^-^-A,^. + -^y(a:,-x,) (1.3) 

j=l i<j 

in the limit of large N. In particular, we are interested in the evolution of factorized initial wave 
functions of the form i/jj^ = ip^^ (here, we use the notation (p'^'^{xi, . . . ,xn) = I\f=i Because 
of the interaction, factorization is not preserved by the time-evolution. However, since collisions are 
very weak, we may still expect that factorization is approximately preserved in the limit of large N . 
In other words, we may expect that, for large N , the solution il)N,t = e~*^"*-0Ar of the Schrodinger 
equation can be approximated (in a sense to be made precise later), by 

^N.t ^ (1.4) 

for a suitable one-particle wave function ipt- Assuming (|1.4p to be correct, it is simple to derive a self- 
consistent equation for the one-particle orbital ft- In fact, (|1.4|) implies that, at time t S R, particles 
are distributed in space, independently of each other, with probability density |(/3tp. This means that 
the potential experienced by a particle at x e can be approximated by the average, mean field, 
potential {V * \(pt\'^){x) and therefore, that 1^94 must satisfy the nonlinear Hartree equation 

idtipt = -Aipt + iV'^\iptf)ipt (1.5) 

with initial data (pt=o = ^■ 

In which sense can we expect the solution of the A'^-particle Schrodinger equation tpN,t to be ap- 
proximated by the factorized wave function on the r.h.s. of (|1.4p ? It turns out that one cannot expect 
convergence in norm (see, however, the recent works [9j [10] where second order corrections to the 
mean-field dynamics are taken into account to obtain a norm approximation of the full dynamics). 
Instead, (|1.4|) has to be understood on the level of the reduced density matrices. Let \ipN,t){'fpN,t\ 
denote the orthogonal projection onto ipN,t- Then, for fc = 1, . . . , we define the fc-particle reduced 
density matrix by taking the partial trace of \ipN.t){'ipN,t\ over the degrees of freedom associated with 
the last N — k particles, that is 

7w,t = Trfc+i^...,Ar \lpN,t){tpN,t\ ■ 

In other words, 7^''^ is defined as the non-negative trace class operator on L'^{M.^'^ ,dxi, . . .dxk) with 
the kernel 

7^''((Xfe;x5,) = J dXk+l ■ . .dxN'^N,t{^k,Xk+l, ■ ■ . , XN)i'N.t{^'k^ ^k+l, ■ ■ -^^n) (1-6) 

where we set = (xi, . . . ,Xk) and, similarly, xj, — {x[, . . . , xj.). From the normalization llV'iv.tll = 1, 
we conclude that Tr^f^^l = 1 for all 1 < fc < iV and all t e R. Observe that, for 1 < A: < A, the 

. (k) 

fc-particle reduced density 7]^ j does not contain the full information about the A-particle system. 

(k) 

Nevertheless, knowledge ofjj^ ^ is sufficient to compute the expectation of A:-particle observables, that 
is of observables of the form O^'^-' (X) l^^"*^) which only act non-trivially on k particles. In fact, 

(OC^) ® l(^-'=)) ^^.,) = Tr7A.,t {o^'^ ® l^^-'^') - Tr7So('=' . 

It turns out that the reduced density matrices are the right quantities to understand (|1.4p . For a large 
class of potentials V, one can show that the reduced density matrices associated with V-'JV,* converge. 
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in the limit of large N, to the reduced density matrices associated with the factorized wave function 
(fif^. In other words, one can show that, for any fixed t e R, 

Tr|7W _|<^,)(^,||^0 (1.7) 

as iV — > oo. Observe that convergence of the one-particle density towards a rank-one projection 
immediately implies convergence of higher order reduced densities as well; for any fixed A; G N and 
t G M, it follows from ^1} that 7^^^^^ ^ \ipt) {ipt\®'' as N — > 00 in the trace norm topology (see Remark[l] 
below) . 

The convergence (|1.7p has first been established by Spohn in [17] for bounded potentials. In [8], 
Erdos and Yau extended the techniques of Spohn to prove (|1.7p for potentials with a Coulomb-type 
singularity V{x) = ±l/|a;| (partial results in this direction were also obtained in 2 ). In \16^, p.7p was 
established again for potentials with Coulomb singularities. In contrast with the previous results, the 
bound obtained in |16| gives an explicit estimate on the rate of the convergence. For factorized initial 
data, it is shown in [16] that 

Tr l^Z-\vt){^t\ <^ (1.8) 

for constants C, K depending only on the initial one-particle wave function ip. The approach used in 
[16] is based on techniques first introduced by Hepp in [12; and then extended by Ginibre and Velo in 
[11] for the study of the related problem of the classical limit of quantum mechanics. More recently, 
bounds of the form (|1.8p on the rate of convergence of the Schrodinger evolution towards the Hartree 
dynamics, were obtained by Knowles and Pickl in [14j, for potential with singularities of the form |a;|~", 
for a < 5/2 (the bound on the rate of convergence obtained in [T3] deteriorates, compared to (|1.8I) . for 
potentials of the form |a:|~", with a > 1). In [3J[T3], the convergence (I1.7P was established for particles 
with a relativistic dispersion (the kinetic energy — Aa;^ is replaced by -y/l — A^^ , for j = l,...,iV) 
and with Coulomb type interaction V{x) = ±A/|x| (this situation is physically interesting because it 
describes systems of gravitating bosons, so called boson stars, and the related phenomenon of stellar 
collapse). In order to describe the dynamics of Bose-Einstein condensates, it is interesting to consider, 
in (|1.3p . two-body potentials which scale with the number of particles N, and tend to a delta- function 
in the limit of large N. In this regime, the many-body quantum dynamics is approximated by the 
Gross-Pitaevskii equation; this problem has been studied in [5j |6j [71 [15] . 

In this paper, we extend the techniques developed in |16) , and we improve the bound (jl.Sp on the 
rate of convergence towards the Hartree dynamics. For interaction potentials with Coulomb type sin- 
gularities and for factorized initial wave functions, we show that the difference between the reduced 
one-particle density associated with the solution of the A^-particle Schrodinger equation and the or- 
thogonal projection onto the solution of the Hartree equation ()1.5p is at most of the order 1/N, for 
any fixed time i e M. The TV-dependence of this bound is expected to be optimal. Note that the same 
bound on the rate of convergence was obtained in 1 4 , for bounded potentials, and, more recently, in 
[1] under the condition that V € L'^(R'^) -t- L°°(R'^), which excludes a Coulomb type singularity. The 
main result of this paper is the following theorem. 

Theorem 1.1. Suppose that the potential V{x) satisfies the operator inequality 

Vixf <Dil-A,) (1.9) 

for some constant D > 0. Let ip G i7"'^(M'^) with \\ip\\2 ~ 1, and let ipt be the solution of the Hartree 
equation 

idm = -A^t + {v*\^t\^)^t (1.10) 

with initial data ipt=o — Let ipN.t — e~^^"*ip^^ and 'y^j^\ be the one-particle reduced density 
associated with tpN.t, o,s defined in (|1.6p . Then, there exist constants C and K, depending only on 
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|(/3||//i and D, such that 



Tr 



7^'l ('^i I (1-11) 



Remark 1. The same techniques used to show il.ll]) can be extended to prove an analogous bound 
for higher order reduced densities; for any fixed fc G N, one can show the existence of constants Ck , Kk 
such that 

Tr\jl,%-\^,)(^r'\<^^- 

Note that, if one is satisfied with a slower rate of convergence for higher order reduced densities, a 
simple argument, outlined in Section 2 of [14 , shows that, for any fc £ N, 



Tr 



7^'l - \^t){^t\''' I < 2^2kTr\^l^l-\^,){^,\\ < 



N 



for constants C,K independent ofk. 



Remark 2. With exactly the same techniques used to prove Theorem \l.l[ one can also consider mean- 
field Hamiltonians with external potential Vext{x) acting on the N particles. The requirements on Vext 
are minimal (the conditions to make sure that the dynamics exists). 

As in [161 [T]. the main challenge to prove Theorem 11.11 consists in controlling the fluctuations around 
the mean-field dynamics. After second quantization, the evolution of these fluctuations is described 
by a two parameter group of unitary transformations U{t; s) (see (j4.5p ). It turns out that the growth 
of the fluctuations can be bounded by comparing first U{t;s) with a simpler approximate dynamics 
U2{t; s) having a quadratic generator £2{t) (see (|4.10|) and (|4.1ip ). A similar approach was already used 
in |16j : to get the optimal bound on the fluctuations, however, we need to consider here, similarly to 
[1], a different approximate dynamics. The problem reduces then to estimating the difference between 
the two evolutions U{t; s) and U2{t; s). While in [T] this difference was bounded using Strichartz-type 
estimates (requiring V £ L^(R'^) + L°°(R^) and therefore excluding Coulomb- type singularity), in the 
present paper we make use of an a-priori bound on the growth of the kinetic energy with respect to 
the approximate dynamics Vl2it; s). It turns out that, for this a-priori bound to be useful, we have to 
introduce a small, A'^- dependent, cutoff ckat in the interaction V] for sufficiently small a^r, we show 
that the error due to the cutoff decays faster than 1/A^ and can therefore be absorbed in the right 
hand side of pTTTI) . 

The paper is organized as follows. First, in Section [2l we show that the many body Schrodinger 
evolution with cutoffed potential remains close to the Schrodinger evolution with full potential V (and, 
similarly, that the Hartree dynamics with regularized interaction remains close to the full Hartree 
dynamics), if the cutoff tends to zero sufficiently fast as N oo. Hence, Theorem 11.11 follows by 
proving the corresponding bound for the difference between the regularized Schrodinger evolution and 
the regularized Hartree dynamics; this crucial bound is stated in Proposition 12.11 In Section [3l we 
define the bosonic Fock space and we recall some of its properties. In Section HI we reformulate 
the convergence problem on the Fock space, and we prove Proposition 12.11 making use of a series of 
estimates (in particular, the a-priori bound for the growth of the kinetic energy, which follows by 
combining Lemma 16.11 and Lemma 16. 2p deferred to Sections [S][51 

2. Regularization of the interaction 
For an arbitrary sequence ajv > 0, we set 

Vix) = sgn(F(x)) • mm{\V{x)\,a],'} (2.1) 
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where sgn{V{x)) denotes the sign of V{x). We also define the regularized Hamiltonian 

N N 

Hn = Y.-A.,,+-Y.^^{x,-x,). (2.2) 

Note that, by definition |V^(a;)| < a^"^ . Moreover, (|1.9p implies the operator inequality 

V'^{x) <D{1- A^). (2.3) 

Instead of proving directly Theorem II. 1[ we show that it is enough to prove the corresponding 
statement for the dynamics generated by the regularized Hamiltonian (|2.2p . if ajv converges to zero 
sufficiently fast. First, we bound the difference between the evolution of the initial iV-particle wave 
function ij^^ w.r.t. and w.r.t. the regularized Hamiltonian H^. 



Lemma 2.1. Let ip_ 



^®JV g^^g ^ ^ H'^{W) with \\ip\\ = 1. Let -tpN.t 



-iH]sit 



N,t 



" ipN- Then there exists a universal constant C > such that 



1pN,t - 1pN,t 



< CNaN \t 



" ^pN and 



(2.4) 



for allN eN,te K. 

Proof. We consider the derivative 



d_ 

dt 



= -2Re — {ipN,t,ipN,t) 

= 2Im {{Hn - HN)ipN,t,i'N,t) 

N 



= — ^Im (^(^V{xi - Xj)- V{xi - Xjij ipN,t,-ipN,t'^ 



i<j 



Observe that the definition (|2.ip of V implies that 

\V~V\<\V\-l{\V\>aJ,^)<\V\''aN. 
Hence, from (|2.5p . we obtain (using also the assumption (|1.9p ) 



d 




2 




di 


IpN.t - i^Nd 




<CN 



{{V{xi - X2) - V{xi - X2)^ ■>pN,t,'ipN,tj 



<CNaN {ipN,u (1 - A.JVa',*)'^' (V'A',t, (1 - A,JVw,t 
Next, we note that, again from ()1.9p . 

A^(^jv,t, (1 - A,J^jv,t) < C{iI;n,u {Hn + N)^l;N,t) 

Similarly, from p.3p . 
Therefore (|2.7p implies that 

The lemma follows after integrating over t 



1/2 



d 




2 




di 






< CNaN 



(2.5) 



(2.6) 



(2.7) 



□ 



As a consequence, we obtain a bound on the difference between the marginal densities associated 
with 'ipN,t and '4'N,t- 
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Corollary 2.1. For any k (z N, let and 'y^j^\ be the k-particle reduced densities associated with 
,/)jv_j = ^^^^ ^ ^-iHNt^(SN^ Suppose aN < in the definition ![KJ\) . Then there 

exists a constant C > 0, independent of k, such that 



Proof. We have 



Tr 



Tr 



(k) ~(fc) 



[k] ~{k) 

i N t - TN t 



sup 

||0<'=)||<1 



< 



N 



(2., 



(2.9) 



where the supremum is taken over all compact operators O'^^^ over L^(M^'^, dxi . . . dxk), with operator 
norm IIOC^)]! < 1. Observe that 

TrO(*^) - In]) = (^n,u (O^'^ <E> 1) ^N,t) - {^n,u {O^'^ ® 1) i^N,t) 

Taking absolute value, we find 

TrOC^) (7^1 - 7^1) I < n^N,t - ^nA (2.10) 
for all observables O'-'"'^ with ||0^'^-'|| < 1. The corollary follows then from Lemma [2.11 □ 

Finally, we estimate the distance between the solutions of the nonlinear Hartree equations with the 
full potential V and with the regularized potential V. 

Lemma 2.2. Let ip G _ff^(]R'^). Let ipt be the solution of the Hartree equation \L10]) and ipt the solution 
of the Hartree equation 

idm = -Alpt + iV '^\(pt\^)(pt (2.11) 
with regularized potential V , with ipt=o — p>t=o = ^- Then 

liV't-^tll <Cawe^l*l. (2.12) 

Therefore 

Tr \\pt){vt\®'' -\^t){^t\®''\<m^t~^t\\<CkaNe''\'\ (2.13) 

for any /c G N. 

Proof. From (|1.9p and (|2.3p it is easy to check that |i<pt||ffi, ll^tllffi < C, for a constant C which only 
depends on We compute 

d 



dt 



2Im (pt, V * \iptY -V * \ipt\ 



= 2Im l^ipt 

+ 2Im {ipt 
where, in the last line, we used the fact that 

Im (^t,(y*(|^t|2-|^,|2))^,^ = 
Using p.6p we find, taking the absolute value, 

^llV't-^tll^ < 2a7v||¥'t||||^J sup / dyV'^{x-y)\Lpt{y)\'^ 



(2.14) 



(2.15) 



+ 2\\Lpt- ^t\\ \Wt\\ sup / dy \V{x - y)\ \ipt{y) - (^t(y)|(|^t(y)| + Wt{y)\) 
< CaN + C \\(pt ~ iptW^ . 



(2.16) 
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In the last inequality we used that, from (|2.3|) . 

dy \V{x ~ y)\ \<ft{y) " ^t{y)\ {\^tiy)\ + \^t{y)\) 



< {^j dy \^t{y) - Uy)?) i^j dyV'ix - y)i\^tiy)\ + ^^'^^^ 

<c\\ipt~f>tuyt\\H^ + \\^t\\H^). 

From (|2.16p we obtain (by Gronwall) 

||<^t-<?t||' < Cajv(e^l*l -1) (2.18) 
which concludes the proof of (|2.12p . To show (|2.13p . we write 

k 

\vt}{^r'' - = E \^t){vr^'-'^ ® i\vt){v^t\ " i^tx^ti) ® \^t){ipr^''-'^ (2.19) 

and we use the fact that 



Tr 



\Vt){'Pt\-\lpt){lpt\ <2\\^t-^t\\. (2.20) 

□ 



As a consequence of Corollarv l2.1l and of Lemma 1^?^ Theorem 11.11 follows from the next proposition, 
which only involves regularized dynamics (regularized A^-particle evolution and regularized Hartree 
dynamics) . 



Proposition 2.1. Let V be as in 12. 1\) with some un > N for some r e N. Let (p e 
7^''^ the one-particle reduced density associated with Q-'^HKt^(g)N solution of the regularized 

Hartree equation 

idt^t = -A^t + iV*\lpt\^)lpt (2.21) 
with initial data ipt=o ~ V'. Then, there exist constants C and K such that 

Tr 7^\) -|^,)(^,| <^. (2.22) 

The proof of this proposition is given in Section U) It makes use of a representation of the problem 
on the bosonic Fock space, which we introduce in the next section. 

3. Fock Space Representation 
The bosonic Fock space over {R^ , dx) is defined as the Hilbert space 

J- = L2(R3, dx)® = " = C ® i2(R3«^ j^^^ 
n>0 n>l 

Here L^(R'^", dxi . . . dxn) denotes the subspace of L^(R'^", dxi . . . dxn) consisting of functions symmet- 
ric with respect to any permutation of the n variables xi, . . . , a;„. In other words, T contains sequences 
if) = {^("^}„>o of n-particle wave functions V'^"^ G Ll{B?'^ , dxi . . . dxn). For V'l, € .F, we define the 
scalar product 



n>l 



A sequence {0, . . . , 0, 0, . . . } describes a state with exactly m particles. We will denote by 

the m-particle sector of J^, which is spanned by vectors of the form {0, . . . , 0, 0, . . . }. The vector 

$7 = {1, 0, 0, . . . } e is known as the vacuum and spans the zero-particle sector J^'^*'^ . 
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The number of particles operator M is defined on the Fock space T by (Afip)^"'^ = n^'^"^. For 
/ e L'^{R^,dx) we define the creation operator a*{f) and the annihilation operator a{f) by 

^''^=1 (3.1) 



(a(/)V')^"^ (a;i,...,a;„) = x/^ y dar /(a;) V<"+i)(a;,a;i, . . . 



For any / e -L^(M^, dec), the operators a*(/) and a(/) are unbounded, densely defined, closed operators. 
The creation operator a* (/) is the adjoint of the annihilation operator a{f) (note that by definition a(/) 
is anti-linear in /), and they satisfy the canonical commutation relations; for any /, 5 € L"^ {M.^ , dx) , 

[a{f),a*{g)] = {f,g)mm), [aif ),a{g)] = [a*{f),a*{9)] = . (3.2) 

For every / e L^{M.^, dx), we introduce the self adjoint operator 

</>(/) = «*(/) + «(/)• (3.3) 

It is interesting to note that the n-particlc product state {0, . . . , 0, f®", 0, . . . } can be produced starting 
from the vacuum fl by applying the creation operator a*{f) for n times. More precisely, we have 

{0,...,0,/®",0,...} = ^^n. (3.4) 



The normalization can be easily checked using the canonical commutation relation. We will be inter- 
ested in the time evolution of these product states. 

We will also make use of operator valued distributions a* and Ux {x gR^), defined so that 

a*(/)= / dxf{x)al 

■'r (3.5) 

«(/) = j dx f{x) 

for every / G L"^ (M.'^ , dx) . The canonical commutation relations take the form 

[ax,a*y] = S{x - y) [a^x,ay] = [al,a*y] = • 

For an operator J acting on the one-particle space L^(R^,dx), we define the second quantization 
rfr(J) of J as the operator on T whose action on the n-particle sector is given by 

n 

(dr(j)v)(") = ^J,^(") 

where Jj = 1^...1^J^1---(S)1 is the operator J acting only on the j-th variable. As an example, 
the number operator is the second quantization of the identity, i.e. Af — dT(l). If the onc-particlc 
operator J has a kernel J{x;y), then the second quantization dT{J) can be written in terms of the 
operator valued distributions dx ) d^ as 



dr( J) = J dxdy J{x; y) a*xay . 



For example, we have 

^ — j dxa^ax ■ 

The following lemma provides useful bounds to control creation and annihilation operators as well 
as operators of the form dr{J) in terms of the number of particle operator Af. 



RATE OF CONVERGENCE TOWARDS HARTREE DYNAMICS 9 

Lemma 3.1. For a > 0, let DiJ\f") ^ {ij; e T : E„>i < 00} denote the domain of the 

operator Af" ■ For any f G L'^(R^,dx) and any G D{M^^'^), we have 

l!a(/)V||< 11/11 IIAA^Vll, 

||a*(/)V|| < 11/11 ||(AA+ 1)1/2^11, (3.6) 

|1^(/)V|| <2||/||||(AA + 1)1/2^11. 

Moreover, for any bounded one-particle operator J on L'^{M?,dx) and for every £ D(Af), we find 

\\dr{J)i:\\<\\J\\\\U^\\. (3.7) 

Proof. These bounds are standard. A proof of p.6p can be found, for example, in Lemma 2.1 of [IB] . 
As for (13.71) . it is enough to observe that 



n>l'Lj = l n>l 

because, clearly, || Ji|| = || J|| for all i = 1, . . . , n. □ 

Given ■0 e J^, we define the one-particle density 7^^^' associated with as the positive trace class 
operator on ^^(M^, dx) with kernel given by 

y) = (^^^^) ■ (3.8) 

By definition, 7^^-* is a positive trace class operator on L'^(R'^,dx) with Tr7^^'' = 1. For an arbitrary 
A^-particle state {0, . . . , 0, ipN, 0, . . . } it is simple to check that p.Sp coincides with the definition (|1.6p 
given in the introduction. 

On T, we define the Hamilton operator Hn by {'Hn'4'Y"'' — ^tv''"'^'""'' ' with 

n 1 ^ 

where V denotes the regularized potential introduces in (j2.ip . Using the distributions ax,a*, Hn can 
be rewritten as 

Hn = J dx\/ xa*^\/ ^a^ + J dxdy V {x - y)a*alayax . (3.9) 

It is clear that, on the A'^-particle sector, the operator coincides with the Hamiltonian Hn defined 
in (j2.2p . To study the dynamics generated by the Hamiltonian FLn on the Fock space, coherent states 
will be useful. 

For / e L'^(R^, dx), we define the Weyl-operator 

W{f) = exp (a*(/) - a{f)) = cxp (^J dx {f{x)al - l{x)a^)^ . (3.10) 
The coherent state with one-particle wave function / is the vector W{f)Q,. Notice that 

W{f)n = e-"/ll'/2 ^ ^^^^f^ = e-"^ll'/2 E {0' • ■ • ' 0' « ■ • ■ > ■ (^-H) 

n>0 ' n>0 ^ ^' 

Hence, if P„ denotes the projection onto the rt-particle sector we have, comparing with p.4p . 

* f -f\n — 11/11^/2 

F„iy(/)1] = e-ll/"'/2^i^^ ^ {0, . . . , 0, /®", 0, . . . } . (3.12) 
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We will use this formula with n — N to write the initial product state (p^^ as the projection onto 
the A^-particle sector of an appropriately chosen coherent state. Eq. ()3.11|) is a consequence of the 
expression 

exp(a*(/) - a(/)) = e-H^H'/^ exp(a*(/)) exp(-a(/)) 

which follows from the fact that [a(/),a*(/)] — commutes with a(/) and a*{f). It implies that 
coherent states are superpositions of states with different number of particles (the number of particles 
is a random variable with Poisson distribution having average . In the next lemma we list some 

key facts about Weyl operators and coherent states. 

Lemma 3.2. Let f,g € L'^{R^,dx). 

i) The Weyl operator satisfy the relations 

Wif)Wig) = Wig)Wif)e^^'''''^f-s^ = Wif + g)e-''"'^f-s^ . 

ii) W{f) is a unitary operator and 

wifY = w{f)-^ = w{-f). 

iii) We have 

W*{f)a,W{f) =a, + fix), and W*{f)a:W{f) = < + J{x) . 

iv) From iii) we see that coherent states are eigenvectors of annihilation operators 

= f{x)i^{f ) ^ a{g)iif) = {g, f)L^^{f) . 

v) The expectation of the number of particles in the coherent state tp{f) is given by that is 

m),mf)) = \\fr- 

Also the variance of the number of particles in ip{f) is given by H/P (the distribution of M is 
Poisson), that is 

m),M^m) - mf)M4'{f)? = ii/f . 

vi) Coherent states are normalized but not orthogonal to each other. In fact 

mf)M9)) = ^ = e-^ll/-^'ll' . 

4. Proof of Proposition 12. II 



Formulating the problem on the Fock space T , the one-particle density associated with the evolution 
g-iffivt^^ of the initial product state ipN = f"^^ has the kernel 

Here we used p.4p . (13.81) and the fact that, on the A^-particle sector, 'Hf^\jr{N) = Hf^. Now, using 
p. 121) . we write 



where we defined 



^ = J,v/9 -P^^(^^^)^ = dN PNW{^/N^)il (4.2) 



- -wi^2 ^ CN'" ■ (4-3) 
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Thus, we obtain 



1 o ^.:W«t„*„ ^-^HNt a*{^) 



N 



-n 



(4.4) 



iv \ 

Following, similarly to [16' , an idea first introduced by Hepp in |12j . we define the unitary evolution 



Hit; s) = W*{VNlpt)e~''^'~'^^"W{VNlps) 



(4.5) 



Lemma 13.21 implies that 



W{t; 0) ayU{t; 0) = W*{VNip)e'^"\ay - V Nlpt{y))e-'^"'W{V Nip) . 



(4.6) 



Moreover, U{t; s) satisfies the Schrodinger equation 



i— W(t; s) = C(t)U(t; s), with W(s; s) = 1 



with the time dependent generator 



dx^ xa*cS/ xttx 



dx (y * \ipt\'^){x) a*aa: 



+ / dxdyV{x ~ y)alayipt{x)ipt{y) + J dxdyVix - y) (^ala*y (pt{x)ipt{y) + a^ay (pt{x] 
^ dxdy V{x - y)al (a*ipt{y) + aj,^t(y)) 



dxdyV{x - y)ala*ayax ■ 



1 

N 



(4.7) 



From (14.41) and (14.61). we obtain 



7fft(2/;^) = ^ /^^l^,M^(ViV^)Z^*(t;0)« + ViV^,(a:))(a, + VN^d.y)p{t;0)n 



(4.8) 



^tiy) (^^^' W^(ViV^)W*(t; 0)a*Z^(t; 0)17 
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Integrating against the kernel J{x; y) of a compact hermitian one-particle operator J (hermiticity 
implies that J{y; x) — J{x; y)), we find 



Tr J ~ |v5t)(v5t|) = y dxdyj{x;y) {j''^\iy;x) - ipt{y) iftix)^ 

N 



m 

N 



n, W{VNip)U*it; Q)dV{J)U{t- Q)VL 

VL, W{^/Nip)U* {t; 0)(t){Jipt)U{t; 0)n 



(4.9) 



where we used the definition p.3p . To bound the second term on the r.h.s. of (|4.9p . we write 

dN /a*{^)^ 



IN 



n,WiVNip)U*(t; 0)(t){Jlpt)U(t; o)n 



In 



N 



n, W{VNip) U;{t; Q)(j}{J^t)U2{t] 0)17 



/N 

dN /a*(y>)^ 



n, WiVNip)Wit; 0) HJipt) {U{t; 0) - U2{t] 0)) Vt 



where we compared the fluctuation dynamics U{t;0) with the dynamics ; 0) defined by the equation 

d_ 



i— W2(i; s) = C2{t)U2{t] s), with ^2(5; s) = 1 
and with the quadratic generator 

L2{t) = JdxVxalVxax+ J dx {V * \(ptf){x)alax 
+ j dxdyV{x - y)alay^t{x)^t{y) 
+ / dxdyV{x-y) (ala*y iptix)ipt{y) + a^rUy (ptix)ift{y) 



(4.10) 



(4.11) 



The existence of the evolution U2{t; s) has been established in [TT]. Taking absolute value in (|4.9p . we 
find 



Trj(7i'i-|^t)(^t| 
dN 



< 



N 



{M + if/^Wit; 0)dr{J)U{t; o)n 



dN 
dN 
dN 



W*(%/jVy>) ° ^'fL n,U*{t; 0) </>( J^t)(Zi2(i; 0) - U{t; 0))n 



(4.12) 
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The first term on the r.h.s. can be bounded, using Lemma 17.11 Proposition 15.11 and Lemma l3. 11 by 



N 



N 



m 



-n 



{Af + 1)^/2 U*{t; 0)dr{J)U{t; o)n 



c 

< — e 

- N 



K\t\ 



Uj\f + 1)^ dr{j)u{t:0)n\ 



< 



^H.ifltl \\{Af + If U{t; 0)n\ 



N 



. C\\J.. 
< — — e 



2^|t| ||(_^+ 1)9/2 <CMgif 1*1 



The second term on the r.h.s. of (|4.12p vanishes. This foliows from Lemma [8. II and because 
PiW*{VN(p)a*{ip)^n = Pi{a*{ip) + VN)^W*{VN(p)n 

The tliird term on the r.h.s. of (|4.12p is bounded, from Proposition 16.21 by 



W*{\fNip) °' ^lL n, {U;{t; O) - U*{t; 0)) <p{Jlpt)U2{t; Q)n 



/N\ 



< 



C||J||e^l*l 
N 



The fourth and last term on the r.h.s. of (|4.12p is bounded, using Lemma [7. II and Lemma ISTT 

dN 

dN 



< 



\N 



< 



{N + lyi'Wit- 0) (t>iJifit)iU2it; 0) - U{t; o))n 



N 



I (A/" + 1) V( J'?t)(W2(t; 0) - U{t; 0))n\ 



Writing {Af+lf(t>{Jipt) = a*{Jipt)iAf + 2f + a{Jipt)Af'^, and using Lemma [511 we find 
dN 



< 



< 



C\\J\\e 



K\t\ 



'N 

C||J||e^l*l 



{Af+lf/^{U2it;0)-Uit;0))n 



N 



where, in the last step, we used Proposition [5TT] 
Summarizing, we showed that 



Tr J 



{inI - l<Pt)(<Pt| 



- N 



for all compact hermitian operators J on L'^{M.^ ,dx). Since the space of compact operators is the dual 
to the trace class operators, and since and \ipt){^t\ are hermitian, we immediately obtain that 



Tr 



^'\-\^t){^t\ 



Tat 



- N 



which concludes the proof of Proposition 12.11 
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5. Bounds on the growth of number of particles 

One of the most important ingredients in the proof of Theorem 11.11 presented in the previous section 
is a bound on the growth of the number of particles with respect to the evolutions U{t; s) and U2(t; s). 

Proposition 5.1. Suppose that lA{t] s) andlA2{t]s) are the unitary evolutions defined in and 
respectively. Then, for every j G N, there exist constants C'j,Kj > such that 

\\{Af+iyU2{t;s)tlj\\ <Qe^^l*'^l||(AA+l)V|| 

and 

\\{M+iyU{t;s)^ < Qe^^^l*-'^l||(AA+l)2J+V|| 

for every ^ Cz J-, t CzR. This implies that, for any j € N, the operators {Af + iyU2{t; s){J\f + 1)^-' and 
{Af + iyU{t; s){Af + 1)^^^^^ extend as bounded operators on the Fock space J- with norm bounded by 

m + iyu2{t; s) {N + i)-^ii < c, e^^i*-^i 

and 

\\{N+iyU{t- s) {N+iy^'^W < Cj e^^'*-^' . 

The proof of this proposition can be found in [16 . More precisely, the bound for the dynamics 
hl{t]s) is given in Proposition 3.3 of 16 . On the other hand, the bound for the dynamics lA2{t]s) 
(which is much simpler), can be obtained using arguments very similar to those of Lemma 3.5 of 
(where a different cutoffed dynamics is studied). 

6. Comparison of Dynamics 

The goal of this section is to estimate the difference between the full fluctuation evolution U{t\ s) 
and the dynamics U2{t] s). 



Proposition 6.1. Suppose that, in the definition V2. 1]) of the regularized potential V , the cutoff a n 
is such that > , for some r G N. Suppose thatU(t]s) and U2{t; s) are the unitary evolutions 
defined in ^.5^ and ^J^Td^, respectively. Then, for any j G there exist constants Cj^r, Kj^r > .such 
that 

m+iy (W(i; s)-Z^2(t;s)) ^11 < ^^^^^4^i^(V',(/C+ AAM2j"+3) + 1)^)1/2, 

V JV 



where JC is the kinetic energy operator 

K. = dr{-A) = J dxV^alV^a^. (6.1) 

Proof. We fix i > and s ~ (all other cases can be treated analogously). Using 

U{t;0)~U2{t;0)^ f dTU{t;T) {£{t) - C2{t))U2{t;0) 
Jo 

we find that 

m+iy {U{t;0)-U2{t;0))ij\\ < f drUAf + iyU{t;T) {£{t) - C2{T))U2{T;0)i;\\ 



< / dr e^(*--) II (AA+l)2J+i£3(r)Z^2(r; 0)^-11 (6.2) 
Jo 

+ /' dr e^(*-") II (AA + ly^+^U U2{t; 0)V|| 
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where we used Proposition 15.11 and we wrote L{t) — >C2(t) — Cz{t) + £4, with 
^3(0 = ^ / dxdyV{x~y)al{al^t{v) + ayVt{y))ax 



1 



(6.3) 



dxdy V(x — y) a*a*ayax ■ 



Using Lemma 16.31 to bound the first term on the r.h.s. of (|6.2p . we find 



< 



c 



-,Kt 



N 



||(AA+l)2j+5/2^|| 



(6.4) 



To bound the second term on the r.h.s. of the last equation we observe that, restricting the operators 
on the n-particle sector J^'^") of the Fock space, 

{n + l)4j"+2 



iV2 



< 



l{n + 1 < iVTJTs) 



< 



N 



N2 



J = l 

1{N +1 > N^-) 



N 



'N 



I.F(") 



where we used the bounds V'^{x) -CI - /\ and \V{x)\ < (see (gH)) and K. is defined in dUT]). In 
Lemma [671] below we show that there exists a constant C > such that, for arbitrary r G M, 

IC < /:2(t) +C(AA+1). 

Hence 

(AA + l)2^-+i d (M + l f=+' < (£2(r) + C(AA + 1)) + ^^-^"^^,-2^"^^ (AA + l 



N 

and therefore 

(W2(r; 0)V, (AA + 1)2J-+1 d iJ\f + 1)2^+1^2 (r; 0)^-) 
1 



iV2c 



< ^ (W2(t; 0)V, (/:2(t) + C(AA + 1))Z^2(t; 0)^) 



(6.5) 



+ T75^(Z^2(r;0)^, 1(AA + 1 > iVW)(AA+ l)4J+6z^2(r;0)^) . 



Using the bound 1(2: > 1) < x"'^'^^+^\ valid for every m € N, Proposition 15.11 and the assumption 
ctN ^ A'^"'^, we can estimate the second term on the r.h.s. of the last equation by 



7V2 



^ {U2 (r; 0)^A,1 (AA + 1 > iV «T5 ) (AA + lf^+'^U2 (r; 0)^) 



< 



< 



^2^^ (W2(r; 0)V^, (AA + l)(4^-+«)("+i)ZY2(r; 0)V) 



J\fm+2-2r 



;^n|(AA+l)(2j+3)(m+l)^||2 
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for appropriate constants C^K (depending on m and j). Fixing m = 2r — 1, we find 

-1^(W2(t; 0)V, 1(A/- + 1 > N^'){N + 1)^^"+«W2(t; 0)V) < -^e^^ W + ■ 

To control the first term on the r.h.s. of (|6.5p we use Lemma 16.21 and Proposition 15.11 together with 
the fact that, at r = 0, 

(V', (£2(0) + C(AA + 1))^) < ilC+M + . 
We conclude that there exist constants C,K > 0, depending on r and j, such that 

C 



Inserting this bound in ()6.4p . we obtain the desired estimate. 



□ 



We also need to bound the difference between the two evolutions U{t;Q) and U2{t\Q) in the third 
term on the r.h.s. of (14.121) . This is the content of the next proposition. 



Proposition 6.2. Suppose that, in the definition \2.1\) of the regularized potential V , the cutoff a n is 
such that aN > N^^ , for some r g N. Let tpt be the solution of the regularized Hartree equation \2.21\] . 
Suppose thatU(t;s) andU2(t;s) are the unitary evolutions defined in |^.5[ j and JQp , respectively. 
Suppose that J is a bounded hermitian operator on L^(R^,(ia;). Then, there exist constants C,K > 
such that 



N 



n, wiVNif) (u*{t; 0) - u^{t; 0)) 0( j^t) ^2(t; o)n 



where dN = y/W./ie-^/^N'^/^) < CN^/* was defined in (J^ 
Proof. We fix t > and we write 



< 



C||J||e^l*l 



dN VN 



n,W{VNip) {U*{t;0)-U;{t;0))<j){Jipt)U2{t;0)n 



ds ( ° n, W{y/Nip)U;{t; s) (£3(5) + C4)U*{s; 0)(t>{Jipt)U2it; 0)n 







where £3(5) and £4 are defined in (|6.3p . Taking the absolute value, we find 



n, w{^/N^) {u* {t- 0) - {t- 0)) 0( j^t) U2 {t- o)f] 



< 



ds 



{N + i)-^'^w*{Vn^) 



{N +lf'''u;{t-s)L:,{s)U*{s-G)(j){J^t)U2{t;Q)^ 



ds 



N 



-n 



{U +lTW{s;Q)(j){J^t)U2{t-,Q)n 



:A + B 



(6.6) 
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where the parameter k > will be fixed later on. To bound the term A we note that, by Lemma l7.1|, 
Lemmal5.ll Lemma 16.31 we have 



U-N Jo 



< 



c 



< 



Nd 
C 



N JO 



y/NdN Jo 



dse^^*""^ \\{J\f+lfU*{s;0)(j){Jiptp2{t;0)n\ 



dse^^lUAf + iy(l>iJipf)U2it;0)n\\ . 



Writing (A^ + 1)^ 4'{J'Pt) = a*{J(pt){J\f + 2)'"' + a{J(pt)Af^, and using Lemma \3A\ we conclude (again 
by Lemma [STT|) that 



A< 



C\\J\\ 



Nd 



N JO 



ds e 



Kt 



(AA+l)"/2^2(i;0)f) 



< 



C\\J\\e 



Kt 



Nd 



(6.7) 



N 



for an appropriate constant K > 0. 

Next, we estimate the term B on the r.h.s. of (|6.6p . On the one hand, we have 

II {M + irW{s; 0)(biJ^tp2it; OM < Ce^' \\ + lf^+'(j){J^tp2{t; 0)fl\\ 

<Ce-^"||J|| iAf+lf^+^/^U2{t;0)n < C|| J|| e-^("+*) 



(6.8) 



where the constants C, K depend on k > 0. To bound the other norm in the term B, we write 



w*iVN^)^^n=J2A., 

y-'^- „i>o 



-n. 



'ml 



Note that, by the unitarity of the Weyl operators (recall that a*{ip)™rt — Vm! {0, . . . , 0, (p'^"^,0, ■ ■ ■}) 



^>l 



Hence 



iJ\f+l)-''C4U2it;s)W*iVNip f*^^ n"^ 



^Nl 



= A,nAi 



m,e>o 

< J2 1-^" 



^M2{i:s){Af +iy Cl{N +l)-''U2{t;s) 



a*{cpy 



-n 



< 



m,e>o 



N^al 



iM+l)-^+'U2{t;s)^n 



E 



^ n,u;it;s)clU2{t;s)^^n 



where < 6 < 1 will be fixed later on. In the regime £ > , we estimated 

1 



4 ^ N^a^ 



■{Af+iy. 



N 



(6.9) 
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In the regime £ < , on the other hand, we used the fact that £4 commutes with J\f, and that 
(TV + 1)^2" < 1. By LemmaEH we have 



< ||(AA + l)-''+2z^2(t;s)(AA+l)«-2| 
1 



(AA + l) 



n 



< Ce^(*-^)/2. 



where the constants C, K depend on n. Therefore, 
1 



E 



2 ^gif(t-.) 
< 



E 



1 Ce^(*-^) 
< 



(£+l)2«-4 - 7V2+A-(2k-6)^ 



To bound the second term on the r.h.s. of (|6.9p . we observe that, on the n-particle sector, 

, / n \ 2 " 



(6.10) 



This imphes that 



•AT 



where K. is the kinetic energy operator defined in ()6.ip . We find 



< 



E 



7V7/4 /l^ \ /^y 



E (^^^,^2*(M)AA^+^^2(M)^f^ 



for arbitrary p > (we use here the fact that \{M > N^^^^) < N p/'^^ J\fP^ for any p > 0). Combining 
Lemma |6. 11 Lemma 16.21 and Lemma l5.ll we conclude that 



E 

£<A* 



< 



E 



^2+p/12q,2 



E 

A ^^jva 



f7, AA' 



4+ 



f7 



< 



7VV4 

iV7/4-25 + jY2+P/12-(5+p)5( 



^2+p/12--(5+p)5^ 



^ ^4+p 



■A 



=A 
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where we used that ||(^t||//i remains uniformly bounded in i e M (by a constant depending only on 
||iy9||/fi). Together with (|6.10p . we obtain, from ()6.9|) . 



^2+5(2k-6)q,2^ ]\[7/4-2S jY2+p/12-(5+p)5Q,2^ 



We fix (5 < f/f2. Moreover, we choose k > so large that N'^+^'^^'^~^^ a% > N^/^, and p > so large 
that A^i+p(i/12-<5)q,2^ > jv3/2_ Then, together with we find that the term B on the r.h.s. of i 

can be bounded by 



B < 



C\\J\ 



IN 



N 



because cIn < CN^/^ (recall that djy = VW./ {e^^^^N'^^^)). Together with dH?]), this completes the 
proof of the proposition. □ 

The next lemma is used to bound the kinetic energy by the generator £2(0 of the dynamics U2{t; s). 
Lemma 6.1. Let £2(0 be as defined in l{4.11\ l and let K, be the kinetic energy operator 



IC ^ J dx Va;a* Vjjfla;. 

Then there exists a constant C > such that the operator inequalities 

-C{N + 1) <£2{t)-IC<C{J\f+l) 

hold true for all i G M. 

Proof. It suffices to show that there exists a constant C such that 

for all ip € F. By definition, we have that 

(-0, {C2[t) - ^ J dx{V * \^t\'^){x){a^tp,a^'ip) + J dxdy V{x - y)^p^{x)iftiy){ay■^p,ax'^p) 

+ 2Re J dxdy V{x - y)!ptix)ipt{y){ayi),al^) . 
The first term on the r.h.s. of the last equation can be estimated by 



(6.11) 



(6.12) 



(6.13) 



(6.14) 



dx(V * |(ptp)(a;)(aa;V',aa;'0) 



The second term on the r.h.s. of (|6.14p can be handled similarly: 



(6.15) 



dxdy V{x - y)(pt{x)(pt{y){ayip, a^ip) 



< / dxdy\Vix-y)my)\'\M\\ 



<lll^l*l^t| 



(6.16) 



dx\\a^i;\\^ < C(V',A/'V')- 
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Finally, the last term on the r.h.s. of (|6.14|) is bounded by 



dxdy V{x - y)ipt{x)Lpt{y){ay'4), altp) 



< J dy\\ay^pr + J dy\^tiyr\\a*iV{y-.)^t)^\\ 
<cfl + sup \\V{y-.)^t\\] {^,{Af+l)ij) 



<C(^,(AA+1)^), 



where we used Lemma |3. II 



(6.17) 



□ 



After controlling the kinetic energy with the expectation of the generator £2(^)1 we have to show 
that this expectation remains bounded in time. This is the content of the next lemma. 

Lemma 6.2. LetU2{t; s) be the evolution defined in ^JTT^, with generator C2{t) . There exist constants 
C and K such that 

\{U2{t;s)4,,C2{t)U2{t;s)i,)\ < Ce^'l*-''! (£2(5) + + l)^') (6.18) 
for all ^ ^ T and all t, s S R. 

Proof. We use the shorthand notation ip2 — l^2{t', s)^- To control {ip2, ■C2{t)'4'2) , we first observe that 

^-{^2,C2m2) = {^2X2^2} (6.19) 

with the time-derivative 

^2(<) = ^ dxdyV{x-y)(^'f^{y)^p^{y) + !p^{y)ift{y)jala^ 

dxdyV{x - y) (^^p^{x)^p^{y) + ipt{x)^tiy)j 0^0^ 



+ 



(6.20) 



+ 2 / dxdyV{x-y)\ipt{x)^Pt{y)ala*y + ^p^{x)ip^{y)a^ay 



Next, we want to control {■4}2^ ^■2{t)'^2) in terms of (^^21 {^2{t) +A/')'02)- There are several contributions 
to the expectation {ip2, C2{t)'ip2) arising from the terms on the r.h.s. of (|6.20p . For example, the 
contribution from the last line of (|6.20p is given by 



(6.21) 



+ 4Im J dxdyV{x~y)ipt{x){V*\(pt\'^){y)ipt{y) {axayip2,il^2) 

since (ft solves the regularized Hartree equation (|2.2ip . Since \\V * |i^tp||oo < ll'^tll^/i !i C, the second 
line on the r.h.s. of the last equation can be bounded by 

dy(y * \^t\^){y) iptiy) {ay^2,a*{V{y - .)(?t)V'2) 



I = 4Re J dxdyV{x - y)ipt{x)ipt{y){ip2,ala*'ip2) 
= -4Im dxdyV{x - y)ipt{x)A^t{y) {axay'il;2,i'2} 



< / dy 



{V * l^tniy) \^t{y)\' \\a*{V{. - y)^t)M' + J dy hy^P; 

< (i + \\v* sup liny - ii(AA+ i)vv2ir 

<C(V'2,(AA+ 1)^-2). 



(6.22) 
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As for the first term on the r.h.s. of (|6.21l) . we write '02 — {'02"^}n>o ^-^d 

Tl>0'^ 

Therefore, introducing the notation x„ = (xi, . . . , x„), 
dxdyV{x - y)lftix)Alftiy) "02) 

= v/(n + l)(n + 2) / dxdydXnV{x - y)(pt{x)A(pt{y)tp^ (x, y, x„) V'2 

n>0 

Integrating by parts, we find 

dxdyV{x - y)ipt{x)Aipf {y) {axayip2,ip2) 

= ^ ViriTWriT2) / dxd?/dx„ VF(a;-2/)(^t(x)V(^t(2/)i/'2"'^^\a;, y,x„)?/'2"''(x„) 

ri>0 

- ^ VT^7+T)(nT2) / dx„y(x- y)^t(x)V(^t(y) Vj,-02"^^V, 2/,x„)0^"^(x„). 

n>0 

In the first term, we integrate by parts once more, but this time w.r.t. the variable x: 
dxdyV(x - y)ipt{x)/S.Lpt{y) {axayip2,ip2) 



ri>0 



n>0 





hl)(n- 


K2) 




hl)(n- 


^2) 




hl)(n- 


^2) 



n>0 

Taking absolute value, and using Cauchy-Schwarz, we find 
j dxdyV{x - y)!pt{x)Aipt{y) {axay^p2,^J2) 

< + J dxdydKnV^ix - y)\^i"+^\x,y,x„ 



n>0 



+ J2{n + 1) fdxdyd^n |V(?t(a;)|2 \W^t{y)\^ |0^"^(x„)P 

n>0 '' 

+ 2^(n + 2) / dxdj/dx„|V,4"^'^(a;,y,x„)|' 

n>0 •' 

+ 2^(n + l) y dxd2/dx„ V^2(a;-y)|^t(x)nv^t(y)nV'^")(x„)|' 
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>2 ^ |,-r.|2|| ^ ll,TC.I|2 



Using the fact that \\V * |(/3t| ||oo < Hv'ill^/i < C for all t e M, and since 



dxdydXnV'^ix - y)\ip!^'^'^\x,y,Xn)\'^ < / dxdydxn ( IV:^')/'^"^''-' (x, y, x„)|^ + (x, y, x„)| 

we conclude that 

dxdyV{x ~ y)ipt{x)Aipt{y) {axay^p2,i'2} < C{ip2, {K. +7V+ l)^p2) 
Together with (|02|) and (jOTj) . this implies that 

|I| <C^(V'2,(/C+AA + 1)V2). 

The contribution from the second line on the r.h.s. of (j6.20p can be bounded analogously. The 
contribution from the first term on the r.h.s. of (j6.20l) is given by 

2Re / dxdyV{x - y)'Pt{y)Vt{y)\\0'x^2\\^ = 2Im / dxdyV{x - y)^t{v)^'Pt{y) ||aa;i/'2|l^ 



, (n+2). 



+ 2Im / dxdyV{x - y)\!pt{y)\'' {V * \vt?){v) Wa^M"" 
A +B. 



The second term can be estimated by 



|B| <2||F*|^t|'llL(^2,AA^2). 
As for the first term, we integrate by parts. Since 

j dxdyV{x~y)\V^t{v)?\\axi:2f 

is clearly a real number, we find 

A = 2Im / dxdy"^V{x - y)^p^{y)V^pt{y) ||axV'2||^ 



Integrating by parts with respect to x, we conclude that 

|A|<4 / dxdy\V{x-y)\\^t{y)\\y^t{y)\\\^xaxM \\a,M 



<4y dxdy\V{x - y)\^ \(pt{y)\^\\'^xaa:H\^ + / dxdy \W'fit{y)\\M2\\ 
< M\V*\Vtf\\oo{^J2,IC^2) + Wt\\m{^J2,Af^2 



<C(^2,(/C+AA)V'2) 

Summarizing, we showed that 

d 



-{lP2,C2{t)^2) 



<C(V'2,(/C+AA+ 1)7^2) 



Together with Lemma [6. II and with Proposition 15. 11 we conclude that 



-(V'2,(/:2(t)+AA+l)^2) 



<C{ij2,{C2{t)+Af+l)^2). 



Hence, the lemma follows from Gronwall inequality. 

Lemma 16.11 and Lemma 16.21 allow us to control the contribution of the term with £4 in 
bound the contribution containing £3 (r) , we use the following lemma. 



□ 

To 
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Lemma 6.3. Suppose C^it) is defined as in i6.3\) . Then there exists a constant C > such that 

\\(j^+iyc,m\<^m+iy+'^'4,\\ 



for all t e R. 
Proof. We compute 

dxdydzdw V{x — y) V{z — w) 

X ('0,a*(a*(^t(y) + ayipt{y))axJ\f'^^ al (a*^ipt{w) +a^^t(w)] a^V) 
2Re / dxdydzdw V{x — y)V{z — w) If tiy)'-Pt{w){'ip, a'^a*ax^J''^^ ala'!i,jaz'i/j) 
+ 2Re / dxdydzdwV{x — y)V{z~w)'ft{y)'fti'^){i^T'^x'^v'^xJ^'^^CLlawaz'ip) 



V 

= 1 +11. 

Using the canonical commutation relations and the formula a^M = {M + l)ax, we find 

I = 2Re J dxdydzdwV{x ~ y)V{z — w)(pt{y)(pt{w){axayazaw 

+ 2Re j dxdydwV{x - y) V{x - w) ^t{y)Vt{w){axayayj{N + ly^^^'^i'.axiN + ly+^/'^ij) 
+ 2Re / dxdydz V{x - y) V{z ~ x) ^t{y)^t{xyaxaya,{M + iy-^'^i>, a,(AA + iy+^/^^) 



= A +B +C. 

(6.25) 

Applying Schwarz inequality, we find 

|A| < j dxdydzdw \\axayazaw{U + iy^'^/'^i)\\'^ 

+ j dxdydzdwV\x~y) V^{z~w)\lft{y)\^\^t{w)\^\\axaz{Af+iy+^/^^^ 

<m+iy+'^'^ijr(i + \\v'*\^t{y)\'\\i 



< C||(AA+l)^+3/2^||2 

where we used that, since < (1 — A), * |^tP||oo < ll'^tllffi ^ C' uniformly in t G R. The second 
term on the r.h.s. of (|6.25p can be bounded by 



|B| < J dxdydw \\axayau,{Af + ly^^^'^i'W^ 

+ J dxdydwV\x-y)V^{x-w)\^t{y)\^\^t{w)\^\\ax{Af+iy+'^^i^\\^ 
< C\\{Af+iy+^^f . 
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Similarly, the third term on the r.h.s. of ()6.25|) is controlled by 

|C| < / dxdydzWa^aya.iAf+iy-^/^tPf 

dxdydzV^ix ~ y) V\z - x) \^t{y)\^\^t{x)\^\\a,{M ^ly^^''^^^ 

where, on the second line, we first integrate over y and we extract the supremum over x of \V'^ * 
Afterwards we integrate over x and extract the supremmxi over z of \{V'^ * |<^(p)(z)| and, 
finally, we integrate over z. 

The second term on the r.h.s. of (I6.24p can be written as 



II = 2Re J dxdydzdwV{x — y)V{z — w) (ptiy)^ti'^){'^xayaz{JV — ly-^jaxawaziJ^ — ^y4') 
+ 2Re / dxdydw V{x — y) V{x — w) (pt{y)'ft{w){axay{Af — lyip, axa-wi-N" — 1)"'V-') ■ 



Hence, we can estimate 

|II| < 2 / dxdydzdwV'^{z-w)\^t{w)\'^\\a^ayaz{M-iyijf 

dxdydw V^{x — w) \ipt{w)\'^ \\axay{J\f — 

< \\v^*\^t\'\\oom+iy+'^'i^f. 

This completes the proof of the lemma. □ 

7. Relation between product states and coherent states 

In this paper we are interested in the evolution of factorized initial data of the form 1^9®-'^ with a 
fixed number of particles N. Since it is more convenient to work with coherent states, we write 

{0, . . . , 0, 0, . . . } = -^%r- = dNPNW{VNv)n 



where the constant d^ — N^^^ takes into account the fact that only a small part of the coherent state 
W{\/N(p)il, lies in the iV-particles sector. Similarly, if we apply the inverse Weyl operator W*{VNip) 
to the factorized state {0, . . . , 0, ip®^ , 0, . . . }, only a small part (of size d^^) of the resulting Fock space 
vector will have a small number of particles. This is the content of the next lemma, whose proof can 
be found in [T]. 

Lemma 7.1. There exists a constant C > such that, for any ip € we have 

C 



< 

dN 



8. A PROPERTY OF THE QUADRATIC EVOLUTION U2{t\ s) 

One of the reasons why we obtain precise error bounds is the observation that, for arbitrary s € R 
and / G L^(IR^, dx), the vectors U2{t; s)a* {f)U2{t; s)Vl and U2{t; s)a{f)U2(t; s)fl are localized in the 
one-particle sector. This fact is proven in the following lemma. 

Lemma 8.1. Suppose that the evolution U2{t;s) is defined as in 114.10^ . Then we have, for any 
f e L2(R3,da;) and any t € R, 

U2{t; Oy^if)U2{t; 0)n = PiU2it; Oy 4>{f)U2{t; m ■ 
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Proof. For any ijj £ J- with \\ip\\ — 1 and ip — 1{M — m)^p, with m 7^ 1, we define the quantity 

= sup 

/eL2(K3) 11/11 

Note that F{0) = 0. Let 



F{t)= sup -]-\{i^,U2it;0ra{f)U2{t;0)n)\+ sup ^ |(V,Z^2(i; 0)*a*(m(i; 0)1^)1 . 
/eL2(K3) 11/11 feL-^iM^} 11/11 



/C = y Va;a* Va;aj, . 
We observe that, for any / € L^(R^,da;), 

Since e~*^* is a unitary operator on ^^(M^, dx), we conclude that 

sup -]-\{^PM2it;0ra{f)U2it■,0)n)\= sup ^ | (V',Zi2(t; 0)*e*'=*a(/)e-*'C*Z^2(i; 0)1)) | 

/GL2(R3) 11/11 /eL2(M3) 11/11 

and similarly if we replace a(/) with a*{f). This implies that 

f^(t) = sup ^ I {i^M2{t- 0)*e*'=*a(/)e-*'CU(i; 0)1)) | 
/eL2(R3) 11/11 

+ sup ^|(V',ZY2(t;0)*e''=*a*(/)e-*'^U(i;0)l))| . 

/GL2(R3) 11/11 

For / e L^(M'^, dx), we compute 

i^MM2{t; 0)*e*'=*a(/)e-*'CU(t; 0)1)) 

= (VAZ^2(i;0)* [e''C*a(/)e-*'CS/:2W-/C]Z^2(t;0)D) 
= {■^,U2{t-Qr [a{ft),C2{t) - K]U2{t;m) 
with /t = e^^^*f. Using the canonical commutation relations, it is simple to check that 
[a{ft),C2it) -IC]^ a{{V * \ij5t\^)ft + {V * fm)^t) + a*{2{V * 'Jm)^t) ■ 
Notice that, under the assumption < C(l — A), we find 

UV*\ipt\^)ft\\ < ||/t||sup / dyVix-y)\^tiy)\' < C\\ft\\\\^t\\lu2 < C\\f\\ 



and 

|2 



* fmm < WtW sup / dyV{x - y)\MyWt{y)\ < q|/t|| < 

X J 

for a constant C, independent of t and /. Therefore, we conclude that 

\{^j,U2{t;0ra{f)U2{t;0)n)\<C f dsF{s) 



11/11 

for every / G L'^{M.^,dx). The same bound can be obtained with a{f) replaced by a*{f). Hence, we 
obtain ^ 

< F{t) <C I dsF{s) . 

Jo 

This bound, together with F{0) — and with the a-priori bound (which follows from Proposition 15. ip 

F{t) < 2\\{Af+l)^/^U2it;0)n\\ < Ce-^l*l 

implies that F{t) = for all teR. □ 
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